Unfrustrated Qudit Chains and their Ground States 
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We investigate chains of d dimensional quantum spins (qudits) on a line with generic nearest 
neighbor interactions without translational invariance. We find the conditions under which these 
systems are not frustrated, i.e. when the ground states are also the common ground states of all the 
local terms in the Hamiltonians. The states of a quantum spin chain are naturally represented in the 
Matrix Product States (MPS) framework. Using imaginary time evolution in the MPS ansatz, we 
numerically investigate the range of parameters in which we expect the ground states to be highly 
entangled and find them hard to approximate using our MPS method. 
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I. INTRODUCTION 



A system with local interactions is not frustrated, if 
the global ground state of the Hamiltonian H = is 
also a ground state of all the local terms each of which 
involves only a few particles. Frustration in a classical or 
quantum system (e.g. a spin glass [lj) is often the reason 
why finding its ground state properties is hard. A locally 
constrained unfrustrated system could still have ground 
states that are hard to find (e.g., 3-SAT where one needs 
to test whether a Boolean formula, made up of 3-literal 
clauses, is satisfiable by an assignment of the Boolean 
variables (§|). 

We choose to investigate chains of (i-dimensional quan- 
tum spins (qudits) with 2-local nearest-neighbor inter- 
actions. Our first result is an analytic derivation of the 
necessary and sufficient conditions for such quantum sys- 
tems to be unfrustrated. Second, we look at their ground 
state properties and find a range of parameters where we 
conjecture that these states are highly entangled and thus 
may be difficult to find computationally. We then cor- 
roborate this by a numerical investigation using a Matrix 
Product State (MPS) method. 

The Matrix Product State description of a quantum 
state has proved to be a very useful tool for the inves- 
tigation of one dimensional quantum spin chains 
A pure state of a system of N interacting d-dimensional 
quantum spins can be written in the computational basis 
as |*) = J2i k =i tp ili2 - tN \h)\i2) ■ ■ ■ \in) with d N param- 
eters For a one-dimensional chain the coefficients 
can be conveniently expressed in a Matrix Product 



State (MPS) form HQ, 
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providing a local description of the system in terms of 
matrices r jfc, [ fe ]. One arrives at this form using a series 
of Schmidt decompositions [8j. The required size of the 
matrices is related to the number x of nonzero Schmidt 
coefficients required for a decomposition of the state into 
two subsystems. In general, \ needs to grow like d N / 2 
for the MPS to be exact. 

Is it possible to capture the essential physics of the sys- 
tem accurately enough with an efficient simulation with 
a much smaller \ ~ poly(-AT), spanning only a small part 
of the full Hilbert space of the system? In our case, the 
qudits are arranged on a ID lattice and only have nearest- 
neighbor interactions. We could thus expect that a re- 
duced space might suffice for our needs. This concept is 
common for the various approaches proposed for efficient 
(tractable on a classical computer) numerical investiga- 
tion of quantum many body systems such as the Den- 
sity Matrix Renormalization Grou p M , Matrix Product 
States @, Tensor Product States |10l| and projected en- 
tangled pair states (PEPSs) [TlT |. 

While the MPS formulation has been shown to work 
very well numerically for most one-dimensional particle 
systems, complexity theory issues seem to show there 
must be exceptions to this rule. Finding the ground-state 
energy of a one-dimensional qudit chain with d = 11 has 
been shown to be as hard as any problem in QMA (l2l . [l3| . 
It is not believed that classical computers can efficiently 
solve problems in QMA. However, to our knowledge until 
now there have not been any concrete examples (except 
at phase transitions) for which MPS methods do not ap- 
pear to work reasonably well. This research was under- 
taken to try to discover natural examples of Hamiltonians 
for which MPS cannot efficiently find or approximate the 
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Figure 1: A qudit chain with nearest neighbor interactions 
H k ,k+i given by ©. The matrices T lk] are a local MPS de- 
scription JTJ of the state. 

ground states. 

The paper is organized as follows. First, in Section 
HT1 we show that the question of non-frustration for qu- 
dit chain Hamiltonians with general nearest-neighbor in- 
teractions can be simplified to only Hamiltonians that 
are sums of projector terms[14]. We then analytically 
show under what conditions zero energy ground states 
for this system exist. Second, in Section HTT1 we show how 
to search for and approximate the ground states numer- 
ically and analyze the efficiency of finding the required 
MPS. We identify an interesting class of unfrustrated qu- 
dit chain Hamiltonians, on which our MPS methods do 
not work well. Led by our numerical work, we conjecture 
that these ground states are highly entangled. Finally, 
we summarize our results and conclude with an outlook 
to further work in Section ITVl 

II. WHEN IS A QUDIT CHAIN 
UNFRUSTRATED? 

We investigate chains of d-dimensional quantum par- 
ticles (qudits) with nearest-neighbor interactions. The 
Hamiltonian of the system, 

JV-l 

H = ^ Hk,k+i (2) 

k=l 

is 2-local (each Hk,k+i acts non-trivially only on two 
neighboring qudits) (Fig. 1). Our goal is to find the nec- 
essary and sufficient conditions for the quantum system 
to be unfrustrated - its ground state is also a common 
ground state of all of the local terms Hk t k+i- The local 
terms can be written as 

H W =E^P^ +1 + J2E^P^ k+1 , (3) 
v 

(k) 

where E is the ground state energy of Hk,k+i and each 
(v) 

Pk k+i i s a projector onto the subspace spanned by the 

(k) 

eigenstates of Hk,k+i with energy Ep . The question of 
existence of a common ground state of all the local terms 
is equivalent to asking the same question for a Hamilto- 
nian whose interaction terms are 

H' Kk +i = Ii,...,fc-i ® Pfe,fc+i ® Ifc+2,...,jv, (4) 

with Pfc,fc+i = J2p=i Pk'l+i projecting onto the excited 
states of each original interaction term Hk.k+i- When 



this modified system is unfrustrated, its ground state en- 
ergy is zero (all the terms are positive semi-definite) . The 
unfrustrated ground state belongs to the intersection of 
the ground state subspaces of each original Hk,k+i and 
is annihilated by all the projector terms. 

We now choose to focus on a class of Hamiltonians for 
which each 

iW=EK*+i>K,*+ii (5) 
P =i 

is a random rank r projector acting on a <i 2 -dimensional 
Hilbert space of two qudits, chosen by picking an or- 
thonormal set of r random vectors (a different set for 
every qudit pair - we are not assuming translational in- 
variance) . 

We now find conditions governing the existence of zero 
energy ground states (from now on, called solutions in 
short). We do so by counting the number of solutions 
possible for a subset of the chain, and then adding an- 
other site and imposing the constraints given by the 
Hamiltonian. 

Suppose we have a set of s n linearly independent solu- 
tions for the first n sites of the chain in the form 

,/.ti,...,i» _ pti,[l]-pi 2 ,[2] pi„_i, [n-l] pi„,[n] f fi \ 

ra n Oil Q1Q2 ' ' * OL n — 2,Oi n — 1 Ot n — l,Oi n \ J 

similar to MPS, with ik = 1, . . . , d and a/. — 1, . . . , Sfe; 
here and below all the repeated indices are summed over. 
The T's satisfy the linear independence conditions 1 

r^'o^a, = 0,Vife,Q!fc_l <S=^ x ak = 0, Va k . 

We now add one more site to the chain, impose the con- 
straint P n , n +i and look for the zero-energy ground states 
for n + 1 sites in the form 

..,i„+i _ ;ii,...,i„_ipi„,[n] pi n+ i,[n+l] ln\ 

The unknown matrix Y 1 ^*^ 1 ^ 1 ^ must satisfy 

(v p I? ?' , n \F*"'N p»n+i.[n+l] _ n (a) 

for all values of a n _i,a n +i and p, with \v„ „ +1 ) vectors 
defined in ([S]) . This results in a system of linear equations 

C pi„+i,[n+l] n (q\ 

with C pan _ uin+1 a n = (w^ n+1 |i„i n +i)r l Q ' l „' [, i, Q „ a matrix 



Note that this is not the standard MPS form, which also requires 
linear independence in the other direction, i.e. 

ya k _ t T l a^\, ak = 0,Vi k ,a k <^=>y ctk _ 1 =0, Va ft _! 

In this case s k would be the Schmidt rank for the partition of 
the qudits into (1, . . . , k) and (k + 1, . . . , n). 
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with dimensions rs n - \ x ds n . If ds n > rs„_i and the ma- 
trix C has rank rs n _i, the conditions (9) are independent 
and we can construct ds n — rs n —i linearly independent 
r^a'I+i" 1 '') corresponding to solutions for the n+ 1 qudit 
chain (see the appendix for further discussion of the rank 
of C) . The freedom we have now is to use only a subset 
of them for constructing solutions. Thus, we obtain the 
formula 

s n +i < ds n - rs n -i, (10) 

valid for all n. The question now is how to choose s n as 
we go along the chain. 

The only constraint on rS, 1 ]'' 1 is linear independence, 
which requires si < d (so = 1 as the first pair have r 
constraints) . If we choose the equality sign in the recur- 
sion relationship above at each step we obtain D n linearly 
independent zero energy states, where 

D n = dD n _i - r-D„_ 2 , (11) 

for all n with Do = 1 and D\ = d. The solution of this 
recursion relation is 

D n = J - J 

f-9 

with f + g = d and fg = r. Hence, 

d [cP d [d 2 

There are three interesting regimes for r and d which 
yield different behaviors of D n (Figure 2): 

1. r > \ gives D n = r 2 with cos6> = ^77. 
D n becomes negative when n + 1 > f and thus 
no zero energy states can be constructed for a long 
chain if r > 

2. r = £ results in D n = (|)™ (n+1), an exponential 
growth in n (except when d = 2, which gives linear 
growth) . 

3. r < 4- implies / > | and / > g so for large n, 

D n ~ /" ^1 — j^j and the number of zero energy 
states grows exponentially. 

Any set of s n that satisfies the inequality (|10[) must have 
s n ^ D n . To show this, we rewrite (|10| as 
s = 1 

si — dsQ = —ui 

s n - ds n -i + rs n -2 = —Un, n > 2 
with m„ ^ 0, n ^ 1. 

These relations can be inverted to give 

S n = D n ~ SLl u l D n-l, Tl ^ 1 

from which s„ ^5 Z?„ follows at once. 




Figure 2: The existence of zero energy ground states for a 
qudit chain with d-dimensional qudits and r projectors per 
pair. We highlight two notable cases: d = 2, r = 1 and d = 
4, r = 4. 

This means that in case 1, it is still not possible to 
construct solutions for a long chain. In cases 2 and 3, we 
can construct sets of states with s„ growing more slowly 
than D n . 

,2 

• For example when d < r < the recursion rela- 
tion (fTU|) can be satisfied also by s n = h n provided 
that h 2 — dh + r < 0. This requires g < h < f, so 
the lower bound on h is the least integer > g. 

• On the other hand, for r < d the simple choice 
s n = 1 also satisfies the recursion (fTU)) . This means 
one can just solve the system from left to right as a 
linear system of equations. This results in a prod- 
uct state solution in the form 

^iii 2 ...i„ = ^[l],ii^[2],i 2 . . . ^[nlin ( 12 ) 

which we can construct by starting with any ip^ 
and finding every from the previous ones. 



III. NUMERICAL INVESTIGATION USING 
MATRIX PRODUCT STATES 

In this Section we numerically search for the ground 
states of our class of random projector Hamiltonians 
We probe the relations obtained in the previous 
Section, and see how well the energy coming from our 
small-x MPS imaginary time evolution converges to zero. 
The numerical technique we use is similar to Vidal's 
[7, 15]. We use imaginary time evolution to bring the 
system from a known state to its ground state: l^grd) = 

limr_s.no ,, c -nJ,^°l I, . In our numerical work we normal- 

\\e WT l*o)ll 

ize the state after every time step [5|. We start from 
a uniform superposition of all the states and Trotterize 
by evolving first the odd pairs of sites and then the even 
pairs. Our experimentation with the parameters for a lin- 
ear chain of length N = 20 is shown in Figures 3 , 4 and 
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Figure 3: (Color online) Ground state energy from imaginary 
time evolution vs. \ f° r different ranks of the Hamiltonian. 
This is a plot for d = 4, and projector ranks of 2, 4, 6. Exact 
description would require \ = d~% = 2 20 . 
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Figure 4: (Color online) This is a plot for d = 5, and projector 
ranks of 4, 6, 8. Exact description would require \ — 5 10 . 



5; all the plots are on semi- log scale and the quantities 
being plotted are dimensionless. 

We see that for r < d the final energy converges to the 
zero energy ground state relatively fast with x "C d N / 2 . 
This can be seen in all the figures by the lowest curves 
(marked by triangles). As can be seen the final energy ob- 
tained from imaginary time evolution tends toward zero 
with a steep slope, indicating that the ground state can 
be approximated efficiently with a small \ m MPS ansatz. 

The r > cP/4 case, marked by squares, is shown by 
the top curves in all the figures. One sees that the final 
energy plateaus relatively fast in all three cases. This 
shows that the numerics have converged to a nonzero 
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Figure 5: (Color online) The case of d = 6, and projector 
ranks: 5, 7, 9, 11. Exact description in general would require 
X = 6 10 . 



value and that increasing \ W1U n °t yield a lower value 
of energy. Therefore, the numerical results suggest that 
there are no ground states with zero energy. 

In the previous section we analytically showed that 
when d < r < cP/4 there are many zero energy ground 
states. However, when we try to numerically find these 
states we see that the final energy converges to zero 
slowly. This is shown in all the Figures by the curves 
marked by circles. Out of these there are the critical 

.2 

cases, where r = These correspond to the curves 
marked by closed circles in Figures 3 and 5. The 
numerical investigation of the case d < r < d 2 /A is 
interesting because it suggests that for large number of 
spins finding the ground state with small Xt tractable 
on a normal computer, is very hard. We interpret this 
as high amount of entanglement among the zero energy 
ground states and leave the analytical proof of this 
statement for a follow up paper. 



IV. CONCLUSIONS 

We have investigated the no-frustration conditions for 
a system of qudits on a line with d states per site and 
random rank r local projector Hamiltonians acting be- 
tween the nearest neighbor sites. We proved that there 

,2 

are no ground states with zero energy for r > and suf- 

.2 

ficiently large N. The system is not frustrated for r < \. 
This second parameter region further splits into two. For 
d < r < —, many entangled zero energy ground states 
exist. On the other hand, for r < d we can also construct 
separable zero-energy ground states (see also Figure |2J). 

We have verified the above numerically, in particular 
we have seen that when d < r < d 2 /4 approximating the 
ground state energy (finding the ground states) is hard 
as the states seem to be highly entangled. Future work 
entails the investigation of the energy gap and the 
amount of entanglement in the system as a function of the 
parameters of the chain. Furthermore, we would like to 
address how far from an eigenstate of the Hamiltonian is 
the wave function after the truncations are made (i.e. as 
a function of \). Finally, we would like to quantify the 
nature of the convergence to the ground state starting 
from an arbitrary state in this framework. 



V. ACKNOWLEDGMENTS 

RM would like to thank Salman Beigi, Michael Artin, 
Sam Gutmann and especially Alan Edelman for fruitful 
discussions. RM and PWS would like to thank the Na- 
tional Science Foundation for the support through grant 
number CCF-0829421. EF and JG were supported in 
part by funds provided by the U.S. Department of En- 
ergy under cooperative research agreement DE-FG02- 
94ER40818, the W. M. Keck Foundation Center for Ex- 



5 



treme Quantum Information Theory, the U.S. Army Re- 
search Laboratory's Army Research Office through grant 
number W911NF-09- 1-0438, the National Science Foun- 
dation through grant number CCF-0829421. DN thanks 
Eddie Farhi's group for their hospitality and gratefully 



acknowledges support from the Slovak Research and De- 
velopment Agency under the contract No. APVV LPP- 
0430-09 and from the European project OP CE QUTE 
ITMS NFP 26240120009. 



[1 
[2 

[3 

[4 

[5 

[6 

[71 
[8] 

[9 
[10 

[11 
[12 

[13 
[14 
[15 



[16 
[17 



K. H. Fischer and J. A. Hertz, Spin Glasses (Cambridge 
University Press, Cambridge, 1993). 

S. Cook, in Proceedings of the Third Annual ACM Sym- 
posium on Theory of Computing (ACM Press, New York, 
1971), pp. 151-158. 

D. Perez-Garcia, F. Verstraete, M. M. Wolf, and J. I. 
Cirac, Quantum Inf. Comput. 7, 401 (2007). 
M. Fannes, B. Nachtergaele, and R. F.Werner, Commun. 
Math. Phys. 144, 443 (1992). 

D. Nagaj, E. Farhi, J. Goldstone, P. Shor, and I. 
Sylvester, Phys. Rev. B 77, 214431 (2008). 

F. Verstraete, V. Murg, and J. I. Cirac, Adv. Phys. 57, 
143 (2008). 

G. Vidal, Phys. Rev. Lett. 91, 147902 (2003). 

M. A. Nielson and I. L. Chuang, Quantum Computation 
and Quantum Information (Cambridge University Press, 
Cambridge, 2000). 

S. R. White, Phys. Rev. Lett. 69, 2863 (1992). 

Z. C. Gu, M. Levin, and X. G. Wen, Phys. Rev. B 78, 

205116 (2008). 

F. Verstraete, M. Wolf, D. Perez-Garcia, and J. I. Cirac, 
Phys. Rev. Lett. 96, 220601 (2006). 

D. Aharonov, D. Gottesman, and J. Kempe, Commun. 

Math. Phys. 287, 41 (2009). 

D. Nagaj, e-print arXiv:0808.2117vl [quant-ph]. 

S. Bravyi, e-print arXiv:quant-ph/0602108vl. 



G. Vidal, Phys. Rev. Lett. 93, 040502 (2004). An alter- 
native yet equivalent notation for the MPS is given in 
[3[. 

M. B. Hastings, Phys. Rev. Lett. 103, 050502 (2009). 
C. Laumann, R. Moessner, A. Scardicchio, and S. L. 
Sondhi, Quantum Inf. Comput. 10, 0001 (2010). 



VI. APPENDIX 



We now address the question of the rank of the matrix 
C in equation (9). This C, which gives the constraints 
on r[™ +1 l depends on and through rW on all 

the \v k k+1 ) with k < n. We would like to say that for 
random \v) and when s k < D k for random choices of 
the s k dimensional subspace, the rank of C is generically 
the maximum rank allowed, min(rs n _i, ds n ). This is not 
obviously true. In particular for the case s k = D k in 
the regime r < to which we restrict ourselves from 
now on, D n grows exponentially in n, while the number 
of parameters in the \v k k+1 ) on which C depends only 
grows linearly. Thus C is far from a generic matrix of its 
size, but we now prove that its rank is indeed rD n —±. 

The argument used by Laumann et al [l7| to prove 
their 'geometrization theorem' also applies to our prob- 



lem. It shows that for a chain of N qudits with ran- 
dom \v p k+1 ), i.e. for a Hamiltonian H as in equations 
(2), (4) and (5) the number of zero-energy states, i.e. 
dim(kcr(iJ)), is with probability one (which is what we 
mean by generic) equal to its minimum value. The 
calculation in section II leading to the recursion rela- 
tion (11) and its solution, shows that in the regime 
r < tf 2 /4 this minimum is > Dn, since if the rank of 
the 7\Dfc-i x dD k matrix C is ever less than rD k —i we 
can choose Sk+i = Dk+i- Hence it is sufficient to find 
a single set of \v k k+1 ) for which dim(ker(_ff)) = Dn to 
prove that Dn is the generic value, i.e. that greater val- 
ues occur with probability zero. This implies that the 
rank of each C is generically rDk-i, since otherwise at 
the first k where C had smaller rank we could construct 
more than D k+ i solutions for a chain of length k+1. 
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V^fc.fc+il^fe+i) = unless 
This can be done for r 

r < 

for d odd are obvious 
Assume that in each Ta^}, ak 
1 to D k . 



the property 
i k s | and i k+1 > |. 
linearly independent \v p ) if 



with 

< 



-p We now assume d is even; the modifications 
We proceed by induction on n. 
with k < n, a k runs from 
From the definition of C (following equation 

+1] is 

This allows us to choose, 



9) and the special choice of 
for i n+1 < "' 
unconstrained for i n +i 







2 and so from equation (9) r^* +1 



< 



2 ' 



for 1 < a 

r; .„ 

with 1 < a 

r 



< %D r . 



1 - 2 



n 

+i,[n+l] 



= 1 when a„ 
< D n , 1 < i n+ i 
= otherwise. 



r = f K - 1) 

< i 



As part of our induction, we assume that for 1 < a n < 



2 D n-l, 
1 a n - 1 ,a n 

with 1 < a 



= 1 when a n = |(a n _i - 1) + i n 



d 



1 < Al-1, 1 < in < 2 

T % an-l,cin = otherwise. 

Now we can show that 



the 



are linearly 



-iMi a » 



independent. 
= for 



rows of 
For if, 
all i n+1 ,a n 



this is true in particular for all i„+i > d/2, 
when it becomes J2 p ypa n -i(K,n+i\ i n'in+i 



<ID 



n-li 



for all 

hi < |; in+i > 21 an d a n-i < D n —\. Since the \v p ) are 
linearly independent, this is only true if ypa„ = for 
all p, a n -\. Hence the rank of C is rD n -\ and a„+i can 
take altogether dD n — rD n ^\ — D n+ i values, which is 
what we wanted to prove. 



